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Abstract 

The axial coupling and the magnetic moment of D*-meson or, more specifically, 
the couplings gD*Dn and <?d*_d 7 , encode the non-perturbative QCD effects describing 
the decays D* — > Dtt and D* — > Dj. We compute these quantities by means of 
lattice QCD with TVf = 2 dynamical quarks, by employing the Wilson ( "clover") 
action. On our finer lattice (a « 0.065 fm) we obtain: g D * Dn + = 20 ± 2, and 
5d*°d°7 = 2.0 ±0.6 GeV -1 . This is the first determination of gr)*°D "/ on t ne lattice. 
We also provide a short phenomenological discussion and the comparison of our result 
with experiment and with the results quoted in the literature. 
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1 Introduction 



Measuring the width of the nearest resonances in the charmed meson spectrum is a challeng- 
ing experimental task. The methods developed by the CLEO-collaboration allowed them to 
make the first measurement of the charged vector meson width, T(D* + ) = 96 ± 22 keV [1]. 
Unfortunately no similar attempt has been made in the experiments at B-factories or at 
CLEO-c. The experimental value for T(D* + ) turned out to be very interesting because it 
provided us with the quantity allowing to check on various theoretical tools that are being 
used to compute the phenomenologically interesting hadronic matrix elements such as those 
needed for the decay constants, form factors, bag parameters and so on. The coupling g c , 
that can be extracted from T(D* + ), as discussed in later sections of this paper, appeared 
to be much larger than predicted by many quark models and by all the techniques of QCD 
sum rules (see ref. [2] for a discussion). The CLEO result was however consistent with the 
old fashioned Adler-Weisberger sum rule combined with phenomenological observations 
made with baryons [2]. 

The first lattice calculation of this coupling in the charmed sector was made in ref. [3] , 
where it was possible to extract this coupling (g c ) in the soft pion limit without running 
into notorious difficulties of dealing with the final state interactions in non-leptonic decays 
on the euclidean lattice. In this paper we provide an update to that result. We use the same 
action and the same methodology as in ref. [3] but the major qualitative difference with 
respect to ref. [3j is that here we get rid of the quenched approximation. The gauge field 
configurations used in this work contain the fluctuations of the Nf — 2 mass-degenerate 
dynamical light quarks. Besides we also implement several technical improvements to make 
the extraction of the form factors from the correlation functions computed on the lattice 
cleaner (double ratios of correlation functions, twisted boundary conditions). On the basis 
of our results we conclude that the lattice evaluation of D* — > decay is consistent with 
experiment, and points towards the upper end of the measured T(D* + ) pQ. 

Beside that quantity, in this paper we report on the first lattice determination of 
the soft photon coupling to the lowest D-meson states, Qd*d^- 111 particular we obtain 
T(D*° — > D°j) = 25 ± 13 keV, in good agreement with experiment, although our error 
bars are still too large to make any stronger statement. The two quantities computed here 
are obviously phenomenologically interesting since the experimental information about 
these decays already exists. Therefore they can be used to check on various quark model 
calculations and to perhaps refine the QCD sum rule analyses. Moreover the quanti- 
ties computed here are particularly interesting to the lattice practitioners. g c -coupling 
is a parameter which appears in the expressions obtained in heavy meson chiral pertur- 
bation theory (HMChPT) that are being used to guide the chiral extrapolations of the 
phenomenologically interesting quantities computed on the lattice, such as the D-meson 
decay constant and its semileptonic decay form factors. Instead, go*D-y is important in 
controlling the impact of the structure dependent term in the radiative leptonic D-meson 
decays in the region of phase space in which the photon is soft. Finally, the radiative 
decay D* — > can be used as a benchmark calculation to compare the lattice approaches 
among themselves and with experiment. The obvious advantage is that in this case one 
goes beyond the hadronic spectrum and computes the hadronic matrix element but there 
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is no need for a CKM parameter to make the precision test of lattice QCD vs. experiment. 

The present paper is organized as follows: we first define the hadronic couplings which 
are the subject of this paper, and relate them to the form factors that are to be computed 
non-perturbatively; in Sec. 3 we explain the strategy to compute the relevant hadronic 
matrix elements on the lattice, the strategy which we then implement in Sec. 4 where 
we also present the results directly accessible from our lattices; in Sec. 5 we provide a 
phenomenological discussion and compare our results with experiment and with various 
model calculations; we finally conclude in Sec. 6. 

2 Definitions 

In this section we define the couplings gD*Dn an d gD*D-yi and relate them to the matrix 
elements that can be computed on the lattice. 

2.1 9d*D-k 

Generically, the coupling of the pseudoscalar (P) and vector (V) mesons to a soft pion is 
defined as 



where q = p — k is the pion momentum and A labels the polarization state of the vector 
meson. Physically this matrix element describes the amplitude of the soft pion emission 
process, V — > Pit, such as K* — > Kn and D* — > Dtt, or the kinematically forbidden, but 
nonetheless very interesting processes, B* — > Bit and uj — > pre. In this paper we will focus 
onto go* Dn, and follow a rather standard procedure, first proposed in ref. jl] and then 
implemented in refs. [U HJ [5]. To describe it in just a few lines we first define the matrix 
element of the axial current = q'j^q (here "g" denotes the light quark, u or d) 



conveniently parameterized in terms of three form factors, Ao^^{q 2 )- We restrain our 
discussion to the charged pion case, so that no anomalous term appears in the axial current. 
At q 2 close to m 2 we can use the reduction formula and write 



(P(k)Tr(q)\V(p,\)) = (e x -q) g VP , 



(1) 



{D{k)\A»\D*{p, A)) = 2m v A (q 2 )^-q» + (m D + m^A^q 2 ) < 



+A 2 (q 2 ) 6X ' q (V + - 




-(D(k)7r(q)\D*(p,\)) = (D(k)\d^\D*(p,X)), 



(3) 



ml — q 



which at q 2 = leads to 



gu*D-K — 
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A (0) . 



(4) 
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This formula is not useful for lattice QCD because the form factor A (q 2 ) is dominated by 
the pion (J p = 0~-state in the t-channel), so that both the q 2 dependence of A (q 2 ) and 
the mass dependence of Aq(0) are extremely steep. However from the fact that no massless 
state can couple to the axial current, we can benefit from the condition that 



2m D *A (0) - (m D * + m D )A 1 (0) - (m D * - m D ) A 2 (0) = 



and eventually arrive at 



9D*Dtt 



m D * + m D 



1 + 



m D * - m D A 2 (0) 
m D * + m D Ai(0) 



(5) 



(6) 



In other words the problem of computing the coupling go*Dn is reduced to the problem of 
computing the form factor A\(0) and the ratio A 2 (0)/Ai(0). We stress again that in the 
above definitions the soft pion is charged go*D-K = 9d*d-k+ • The coupling to the neutral 
pion is related to the one we compute here via isospin, i.e., g 2 D , Dn + = 2>9d*d-k°- 



2.2 go*Dj 

The coupling of the pseudoscalar (P) and vector (V) mesons to a photon is defined as 

( 7 (<Z, r)y)P(k)\V(p, e x )) = e v^x u p a k? g VPl , (7) 

which we will accede by computing the matrix element of the electromagnetic current, 
namely 



{P(k)\JT\V(p,e x )) 



6 SiivaB £\P kP ■ 

my + VTip 



where Jf™ 

A 1 



QqQi^iQ + QqQl^Qi an d the matrix element is expressed in a rather standard 
way, i.e. in terms of the form factor V(q 2 ). In the above expressions e = y/4ira cm , and 
Qq,9 is the charge of the heavy quark and the light antiquark (or vice versa). It is very 
important to keep track of the relative sign difference between the two terms in J° m . In 
the case of D-mesons we have 



: F d (0) , with F d {0) = lv qq (0) 

mo + mo* 3 

F u (0) , with F u (0) = lv™(0) 

mo + mo* o 



1 
2 

1 + 



V cc (0) 
+ W«(0) 
V cc (0) 
Vii (0) 



(9) 



We evidently separated the 'charm' (cry^c) from the 'up/down' {q^y^q) contributions to the 
electromagnetic current the matrix elements of which will be computed separately. We 
should emphasize that we do not consider the isospin violating effects and we consistently 
take m,, 



m d 



m 



Notice that, contrary to the pionic coupling (go*Dir), the coupling to 
the soft photon (go*Dj) is dimension-full, and it is a measure of the magnetic moment of 
the D*-meson. To make the lattice computation more straightforward our target will be 
to compute the dimensionless form factor F u d (0). 
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3 Correlation functions computed on the lattice 



In this section we list the correlation functions that are to be computed on the lattice, 
define the convenient double ratios of correlation functions which lead to the desired form 
factors. We then discuss the method that helps exploring the kinematical configurations 
in which when q ^ 0. 

3.1 Axial form factors 

In order to reach the matrix element (|2J) we need to compute the following three-point 
correlation functions 

where V M = crf^q, and V = crf^q are the interpolating operators for the vector and the 
pseudoscalar meson respectively, while Vj* = q'j^q and = qj^q. The sources are 
fixed at t — and ts = 24, and the matrix element is extracted from the usual time 
dependence of the correlation functions. In recent years it became clear that a more 
reliable information can be extracted if one cancels the exponential time dependencies in 
the correlation functions by combining them into suitable double ratios [6j [3, Ej . This is 
why in eq. (TTUj) we also defined the elastic correlation functions, E^- 1 ' 2 '. In the case in 
which both D- and D*-mesons are at rest the only contributing form factor is A^q^^), 
with g^ ax = {m D * — m D ) 2 . In that case, the convenient double ratio to consider is 

(D(0)\A \D*(0))(D*(6)\A \D(0)) 
(D(fy\V \D(fy)(D*@)\V \D*(d)) 

[(m D . +m D )A 1 ( y ql^)] 2 

Since in eq. (jSj) we need the form factors at q 2 = 0, an extra step is needed. As it is 
well known, on the periodic cubic lattice of size L, the smallest momentum that can be 
given to a particle is g m i n = (2n/L), which on the currently accessible lattices would be far 
too large, and would eventually push g 2 's to relatively large negative values. It has been 
shown in refs. j9] that by imposing the twisted boundary conditions on one of the quark 
propagators in the correlation function, one can explore the arbitrary small momenta of 
the hadron associated to that "twisted" quark. A simple recipe to implement this idea in 



= ^(V M (0,0)^(f,t)pt(^t 5 ) e -^) , 
= XJ(P(0, 0)A fl (x,t)Vt(y,t s )e- i ^) , 
= J^(V(6,0)V q (x,t)V\y,t s )e-^) , 

x,y 
i,x,y 



= c&\d-,t)cs\d-,t) 
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practical calculations consists in rephasing the gauge links as 

UM -> U fa) = e^U^x), where 0„ = (0, 9) 



and then inverting the propagator on the lattice with periodic boundary conditions, i.e. 
by inverting the Wilson-Dirac operator on such a rephased gauge field configuration, 
S q (x, 0; U e ) = (q(x)q(0)). Finally the twisted propagator is obtained as 



(12) 



The net effect is that the resulting ground state extracted at large time separations between 
the sources in two-point correlation functions, e.g. 



^<Tr[sf(0,x; U) l5 S q (x, 0; U) l5 \ ) 

X 

satisfies the following dispersion relation 0, [10 

El 



2E, 



■|(0|975?|t> 



(13) 



V- 



(14) 



Therefore we can give as small a momentum to a daughter pseudoscalar meson in eqs. (|2|8|) 
as we want. An important discussion and description of the same physics by means of a low- 
energy effective theory in ref. [H] assert that adding a small twisted angle to the boundary 
condition on the valence quark, but not on the sea quarks, has a negligible impact on the 
low energy QCD observables, such as fx,n and m^, and we shall neglected it in what 
follows. For further reference, when using the twisted quark propagator, we take vector 9 
to be of the form 9 = (9q, 9$, 9q). We will also always keep the decaying vector meson at 
rest. More specifically, we compute the correlation functions ( TTUl) as follows 



E^(q-t) 



(J2 Tr [S c (y, 0; U)^S q (0, x; U) lul5 sfa y- C/) 75 ] ) 
(J2 Tr [S c (y, 0; U) l5 S°(0, x; U) llllh S q {x, y; U) lu ] ) 
(J2 Tr [S c (y, 0; U) l5 S°{0, x; U) lo sfa y- t/) T5 ] ) , 

x,y 

o Tr ^ c ^' 0; U )liS q (0, x; U)^ S q (x, y; U)^}) , 



(15) 



i,x,y 



where q = 9/L, and (...) denotes the average over gauge field configurations, U. To extract 
the form factor at \q\ ^ 0, we first note that for large time separations one has 



i=l 
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Figure 1: The Feynman diagrams used to compute the ratios R\ (t) and R% (t) in cqs ^^) and 
respectively. The full dot symbols are used to denote the transition operator, the matrix element of which 
we are interested in. The symbol "6" with a given quark line indicates that the twisted boundary conditions 
were imposed in the propagator inversion of that particular quarJc. The single- and double-line distinguish 
between the light and the heavy quark propagator. 



PWm,-** x (run + m D .) M f) x «1^M , ( 16 ) 

and therefore the convenient double ratio to extract Ai(q 2 ) looks very similar to Ro(t) in 
eq. ffTUl) , and reads 

R (t) = C^m^fm ^ [(m D .+m g )A(g 2 )] 2 

With our D*-meson at rest the point q 2 = is reached when the twisting angle 

6' = J_ m2 D*-< . (18) 

In practice, we give a few values of 9q around 9' so that we can interpolate the form factors 
to q 2 = 0. As for the correcting factor ^(g 2 ) / A\(q 2 ) we extract it from the following ratio 
of the correlation functions 



, E D - m D * ri(2) t 

2qim D , A 2 (r 



m) effm "k+s 1 *)' (19) 

where, with the definition (TTTj) . q { = 9 /L. This concludes our strategy for computing the 
axial coupling, gD*Dn, on the lattice. 
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3.2 Vector form factors 



Concerning the coupling to a soft photon, we should first emphasize that our vector meson 
is at rest so that the matrix element (jHJ) simply becomes 

(D(k)\J° cm \D*(0,e x ))=0, and (D(k)\ J cm \D*(0, e A )) = (e A x k) ^ ™ D * V{q 2 ) .(20) 

m D * + m D 

The correlation functions needed in this case are 

C\] qq \q; t) = Tr [S e (y, 0; 17)^3,(0, x; U) l3 S S q {x, y; U) l5 ] ) , 

x,y 

Ci^m = -<^Tr[5 c (a;,0;t07*5 g (0,y;r076Sf(y,a:;t07i]>, 
C% qq) (q;t) = (^Tr^ c (y,0;C/)75^(0,x;C/)7^(x,y;^) 7i ]>. (21) 



Notice in particular the sign difference of "eye"- with respect to the "gyV'-part of the elec- 
tromagnetic current, which comes from the Wick contraction. Physically that accounts for 
the spin difference between the D and D* states. Finally, that minus-sign compensates the 
relative sign difference between the electric charge of the light and heavy quark/antiquark 
in the heavy-light meson. That sign-flip has been properly taken into account when writing 
the expressions (Q. Finally from the double ratio 



>(2<?<Z) 



R\{t) 



EW(q;t)EW(0;t) 



3 E D (m D * + m D 



1 2 ' 



(22) 



we extract the desired form factor. In the above expression we use the symmetry of the 
problem and average over the equivalent indices, e.g. 



c% qq) m 



C[l qq \q ] t)+Ci qq \q ] t) + Cif q \q ] t) 
C^ qq \q;t)-Cg qq \q-,t)-Cil qq \q;t) 



(23) 



Notice yet another "— " sign in eq. (1221) which is meant to compensate the sign difference 
between the two matrix elements accessed from (q; t) (photon emission, D* — > D7) 

and from q9 \q;t) (photon absorption, 7Z) — > D*). Finally, the ratio of the two form 
factors is computed through 



R v 2 {t) 



Q] 9q \q;t) VM{q*) 



(24) 



where Clj CC \q;t) is of the same form as C^"\q;t) after replacing q <-> c in eq. (121]) . The 
Feynman diagrams of the double ratios leading to R\{t) and R\{t) are shown in fig. [T^i 
and [lb respectively. 
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c sw m 




z v nans] 


z A USE! 


b v m 


b A pt - m 


5.29 
5.40 


1.919 
1.823 


-0.034 
-0.032 


0.743 
0.757 


0.772(4) 
0.783(4) 


1.91 
1.79 


1.31 
1.30 



Tabic 1: Rcnormalization and improvement constants used in this work. Apart from cy and b A , which 
are estimated by using the 1-loop (boosted) perturbation theory, all the constants are determined non- 
perturbativcly. 



3.3 Improvement and renormalization 

The gauge field configurations that are being used in this work are obtained with the 
Wilson plaquette gauge action, while the effects of dynamical quarks are incorporated by 
using the 0(a) non-perturbatively improved Wilson fermions, with "a" being the lattice 
spacing. We consider only the fully unquenched situations in which the light valence and 
the sea quarks are mass degenerate. To fully respect the 0(a) improvement we also need 
to improve the operators. The subtraction of 0(a) effects needs to be done at the level of 
bare lattice operators, and through the mass correction to the renormalization constants. 
In short, 



4 mpr -(x) = Z A (g 2 )(l + b A (g 2 )(am q )) [A,(x) + c A d,P(x 

V^-(x) = Z v (g 2 )(l + b v (g 2 )(am q )) [V,(x) + c v d u T, v {x)] , (25) 

where T^ u (x) = iq(x)a tll/ q(x), and o^ lv = (i/2)[y li ,j u \. When sandwiched between the 
external V and P states, or more specifically D* and D, the piece proportional to ca{Qq) 
does not modify the form factors Ai^iq 2 )- It only changes Ao(q 2 ) as 

^ r V)=(l-c^o 2 )^) , (26) 

where m^ WI stands for the bare quark mass obtained on the lattice via the axial Ward 
identity (AWI). In other words, our determination of the coupling go*Dn does not feel the 
effect of improvement of the bare axial current. Instead, the vector form factor does get 
improved as 

V i ^iq 2 ) = V(q 2 ) + c v (g 2 )(m D + m D *) T^q 2 ) , (27) 
where we used the standard parameterization of the tensor density matrix element, 

(D(k)\TnD*(p,e x )) = {e^ [L + k p - ^"^ J Uq 2 ) + ^ q,Uq 2 ) 
+^e^k aP ' p (r 2 (q 2 ) - Tl (q 2 ) + ^ ml Ts^)) } e X a (p) . (28) 

We implemented this improvement in our calculation and the results presented in the 
next section are improved and renormalized. The values of the improvement (cy(<7g)) and 
renormalization constants (Zy^S'o), bv,A(9o)) used in this work, with the appropriate list 
of references, are listed in table [U 
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4 Numerical results 



4.1 Lattice details 

The correlation functions needed to complete this work are computed on the gauge field 
configurations produced by the QCDSF collaboration [19]. They were obtained by using the 
Wilson plaquette gauge field action, and the improved Wilson fermions with Nf = 2 mass- 
degenerate dynamical light quarks at /3 = 5.29 (a ~ 0.075 fm) and (3 = 5.40 (a ~ 0.065 fm) 
on the lattices of volume 24 3 x 48. From the ensemble of available configurations we chose 
the ones corresponding to the lighter sea quarks, the hopping parameters of which are 
listed in table [2] The number of configurations used in this study is also given in the same 
table. We chose to use the configurations that are separated by 20 molecular dynamics 



p 


^sea ^val. 


# conf. 




m D 




5.29 


«i = 0.1355 


60 


0.3271(23) 


0.789(3) 


0.854(4) 




k 2 = 0.1359 


80 


0.2451(25) 


0.760(5) 


0.805(6) 




k 3 = 0.1362 


100 


0.1549(24) 


0.731(6) 


0.783(9) 


5.40 


«i = 0.1356 


130 


0.3124(18) 


0.725(3) 


0.771(3) 




k 2 = 0.1361 


120 


0.2166(23) 


0.691(4) 


0.734(5) 




k 3 = 0.13625 


160 


0.1843(26) 


0.681(2) 


0.721(4) 



Table 2: Masses of the light-light and heavy-light mesons in which the light valence quark has the same 
mass as the sea quark. All results are given in lattice units and the charm quark hopping parameter is fixed 
t° K Ihami 25 ' = 0-125 and ^|£~ 5 m 40 ' = 0.126. The gauge field configurations are separated by 20 unit-length 
HMC trajectories. 

trajectories of the unit length in the hybrid Monte-Carlo algorithm (HMC). In computing 
the quark propagators we attempted using the so called "color dilution" technique [21] 
but since we did not observe any appreciable effect in improving the statistical quality of 
our correlation functions, we returned to the standard BiCG-Stab algorithm [22] • With 
the computed propagators we then calculated the 2- and 3-point correlation functions in 
a standard way. We first recomputed the masses of light quarks and of the pseudoscalar 
mesons consisting of the light valence quark and antiquark the mass of which is equal to 
that of the sea quark. 

We tried to fix the charm quark mass in several ways, of which we describe one. Since 
the charm quark is not propagating in the sea, its phenomenology is quenched and its 
hopping parameter can be fixed by tuning its value to make the ratio of the light-strange 
pseudoscalar and the strange-charmed mesons equal to its physical value, m ria Jmr, s = 0.35, 
where m Vss = 2m? K — m\ = 0.684 GeV. On the very same configurations the strange 
quark mass has been computed in ref. |18| . quoting rom^ ls (2 GeV) = 0.2545(47) at (3 — 
5.29, and r m^ s (2 GeV) = 0.2671(48) at (3 = 5.40. By using these values and the numerical 
constants provided in ref. [TS], as well as the force parameter [2D], r /a = 6.25(10)^=5.29, 
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7.39 (26) / 3 =5 .40; we extract Strange ~ 0.1355 and ^range^ ps 0.1359. After having fixed one 
valence quark to ^strange and varying K c harm, we were able to reproduce m Va Jmr, s = 0.35, 
for ^chamf 5 ' = 0.125, and ^jf^f = 0.126, after linearly extrapolating the sea quark mass 
to zero. By adopting other strategies to fix the charm hopping parameter, the resulting 
ftcharm differs by at most 5 x 10~ 4 . We checked, however, that the results of this paper are 
insensitive to that variation. In table[2]we list the pseudoscalar (D) and vector (D*) meson 
masses, for which the light quark hopping parameter K q = n sea and the heavy quark is fixed 
to ft c harm- Notice the labels ^1,2,3 in table [21 which will be used in the next subsections 
when reporting our results for the form factors. Every time the light valence quark appears 
in a correlation function, its mass is always kept equal to the sea quark mass. 



4.2 Form factors 

In this subsection we present our main numerical results. We first computed the ratios (TTU|) 
which exhibit the long plateaus and lead to y4i(g^ ax ), with q^^ = (mo* — mo) 2 - Since 
<?max 1S very small and since the shape of the axial form factor is expected to be driven 
by a rather heavy ai-state, the result for Ai(0) will differ only slightly with respect to 
^i(<?max)- 1 To control the interpolation to q 2 = we evaluate the ratios (fT9l by chosing 
the twisting angle 9 G {0.5, 1.0, 1.5}, which are sufficient to make q 2 straddle around zero 
for all of our ^1,2,3- 

From our data we observe that the quality of the signal for the ratio R\(t) remains 
good when 6q 7^ and small, whereas the signal for R 2 (t) is bad. The illustration of the 
ratios R\^{t) is provided in fig.[2J Whenever the contribution of the form factor A 2 (q 2 ) to a 
correlation function is significant, those correlation functions are numerically much smaller 
than those in which A 2 does not contribute, e.g. jC^ (q;t)\ <C (q;t)\. As a result 
the extracted information on A 2 {q 2 ) suffers from large statistical errors. We attempted 
several other options to extract A 2 (q 2 ) / Ai(q 2 ) but none appeared to be better than the 
one based on using eq. (1191) . the results of which we present here. Luckily, however, the 
ratio A 2 (q 2 ) / Ai(q 2 ) in eq. comes only as a small correction to the dominant Ai(q 2 )- 
contribution, i.e. it is suppressed by {m D , — m D )/{m D , +m D ). In table [3] we list our 
results for the form factor A\{q 2 ) as obtained from the fit to Ro t i(t) to a constant between 
t G [10, 15]. The results for A 2 (q 2 ) / A\(q 2 ) are combined into the quantity 



2 2n 2y/m D m D , 2 

9c{q ,mj = A x {q ) 

m D + m D * 



1 mp* - m D A 2 (q 2 ) 
m D * + m D Ai(q 2 ) 



(29) 



that will eventually lead us to the g c -coupling, defined as, 



9d*d-k = 7 9c ■ (30) 



1 From the light-light two point correlation functions, CAiAiif) = Sz^iW^M 2 -) we extract the masses 
of ai-state and from our data we have: 

{m a Jm v ) = {1.9(2), 2.5(3), 3.9(3)^=5.29, and {m ai /m n ) = {2.0(0), 2.6(1), 2.9(l)} /3=5 .40- 
On our lattices, m ai is about 10 times larger than g£ 
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Figure 2: In the upper plots we illustrate the time dependence of the ratios Ro(t), Ri(t) and R^t), 
defined in egs. lfl0\17\19\) . leading to the dominant form factor A\(q 2 ) and to A2(q 2 )/Ai(q 2 ). In the lower 
plots the similar illustration is provided for the ratios R\ (t) and RX(t), defined in eqs.|22|) and |24|) . 
and leading to the form factor V qq (q 2 ) and V cc (q 2 ) /V qq (q 2 ) , respectively. The fit on the plateau interval 
t e [10, 15] is also shown, for both 9 = 0.5 and 9 = 1.0, corresponding to about 110 MeV and 220 MeV 
respectively. 



Concerning the momentum injections corresponding to 6q G {0.5, 1.0, 1.5}, by chosing r$ = 
0.467 fm, we have \q |^ =5 . 29 G {96(2), 191(3), 287(5)} MeV, and \q |^ =5 .4o G {113(4), 226(8), 
339(12)} MeV. These values would be 7% smaller if we made another (also standard) a 
choice, namely r = 0.5 fm. The results for all our data-points (both j3, all K sea , and all 8 ) 
are collected in table [3j Notice that at /3 — 5.40 we did not compute the form factors with 
9q = 1.5, as it was clear that for that value q 2 <C and would not help us interpolating to 
q 2 = 0. 

We now turn to the results for the vector form factor. In this case the form factor 
is not accessible when both mesons are at rest. When # ^ we obtain the results for 
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/3 = 5.29 


p = 5.40 




{r q) 2 


AM 2 ) 


g c (q 2 ,ml) 


(r q) 2 


Ai(q 2 ) 


g c (q 2 ,ml) 




0.17(2) 


1.40(11) 


1.45(11) 


0.12(1) 


1.11(4) 


1.15(4) 




0.09(3) 


1.06(10) 


1.09(11) 


0.11(2) 


0.98(7) 


1.01(7) 


K-3 


0.10(3) 


0.93(14) 


0.95(15) 


0.09(1) 


0.90(3) 


0.93(3) 


Kl 


0.11(2) 


1.37(11) 


1.42(19) 


0.04(1) 


1.09(4) 


1.10(7) 


K-2 


0.03(3) 


1.04(10) 


0.98(17) 


0.03(2) 


0.96(7) 


1.04(2) 


^3 


0.04(3) 


0.89(13) 


0.88(26) 


0.02(1) 


0.88(3) 


0.91(13) 


Kl 


-0.05(2) 


1.29(10) 


1.34(11) 


-0.19(1) 


1.03(4) 


1.05(4) 


K-2 


-0.13(2) 


0.96(10) 


0.99(13) 


-0.20(2) 


0.87(7) 


0.92(8) 


^3 


-0.12(3) 


0.85(18) 


0.83(17) 


-0.21(1) 


0.79(3) 


0.80(3) 


Kl 


-0.33(2) 


1.16(10) 


1.21(10) 


-0.56(1) 






«2 


-0.40(2) 


0.85(9) 


0.88(11) 


-0.57(1) 






^3 


-0.39(3) 


0.66(14) 


0.64(12) 


-0.58(1) 













V m {q 2 ) 


V cc (q z ) 
Vii{q 2 ) 


(r qf 


V qq {q 2 ) 


V cc {q') 
Vii{q 2 ) 


O = 0.5 


Kl 


0.11(2) 


3.73 ± 1.63 


0.03(39) 


0.04(1) 


5.03 ± 0.59 


0.50(9) 




K2 


0.03(3) 


2.80 ± 1.39 


0.89(48) 


0.03(2) 


4.83 ± 0.88 


0.42(9) 




^3 


0.04(3) 


5.69 ± 3.34 


0.32(45) 


0.02(1) 


4.35 ± 1.12 


0.42(17) 


e = i.o 


Ki 


-0.05(2) 


4.82 ± 0.89 


0.34(8) 


-0.19(1) 


4.73 ± 0.34 


0.50(5) 




K2 


-0.13(2) 


4.38 ± 0.74 


0.51(10) 


-0.20(2) 


4.48 ± 0.59 


0.46(5) 




K-3 


-0.12(3) 


6.07 ± 2.02 


0.35(20) 


-0.21(1) 


4.15 ±0.64 


0.42(8) 


O = 1.5 


Ki 


-0.33(2) 


4.76 ± 0.68 


0.41(6) 


-0.56(1) 








K-2 


-0.40(2) 


4.45 ± 0.55 


0.47(7) 


-0.57(1) 








K-3 


-0.39(3) 


4.70 ± 1.00 


0.38(10) 


-0.58(1) 







Table 3: The form factors Ai(q 2 ) and g c (q 2 ,rn 2 ) [c.f. eq. ((29|) ] relevant to the axial coupling, and the 
form factors V qq (q 2 ) and V cc (q 2 ) /V qq (q 2 ) relevant to the magnetic moment of our D*-meson. 
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(3 = 5.29 


(3 = 5.40 










Ki 


K2 




Vii (o) 


4.56 ± 1.12 


3.14 ± 1.43 


5.88 ±2.88 


4.98 ±0.54 


4.78 ±0.82 


4.34 ± 1.08 


V cc {0)/V qq {0) 


0.26 ±0.13 


0.76 ±0.36 


0.33 ±0.36 


0.50 ±0.08 


0.43 ±0.08 


0.42 ± 0.16 


K{0,ml) 


3.62 ± 1.17 


3.71 ± 1.57 


5.19 ±2.14 


4.98 ±0.45 


4.55 ±0.81 


4.10 ± 0.89 




-0.95 ±0.63 


-0.56 ±0.77 


-0.69 ± 1.55 


-0.00 ±0.26 


-0.23 ±0.27 


-0.24 ± 0.52 


9c(0,ml) 


1.37 ±0.11 


1.02 ± 0.12 


0.89 ±0.14 


1.10 ±0.04 


0.99 ±0.07 


0.90 ±0.06 



Table 4: Results of the linear interpolation of our results from table [3| to q 2 = 0, for each of the light 
quarks that we have at our disposal (Ki t 2.z)- We also show the particular combinations F U) d(0, to 2 ) which 
are needed to arrive at the magnetic moments, i.e. to the couplings .9_d*+d+ 7 and g D *o D o 7 . 



V qq (q 2 ) and V cc (q 2 )/V qq (q 2 ) from the fit of the ratios R\{t) and R%(t) to a constant, on 
the plateau which we chose to be tfi t € [10, 15], as before. The quality of the signal is 
shown in fig. [2, and the complete list of results is provided in table [31 From that table 
we see that, contrary to what one would naively expect, the charmed contribution to the 
electromagnetic current is not negligible when compared to the light quark piece. In other 
words around 30% of the magnetic moment of the D*°-meson comes from charm. That is 
an indication that indeed for some quantities, such as the case at hand, the leading term 
in heavy quark expansion is a bad approximation of the charmed hadrons properties. 2 

4.3 Towards the physical results 

To reach the physical results, we need to make two important steps. First we need to 
interpolate to q 2 = 0, which we can do either linearly (or quadratically) , or by using some 
physical assumption, such as the nearest pole dominance of the particle exchanged in the 
t-channel, namely ai-meson to g c {q 2 , m n)> an d the p-meson dominance to the V qq (q 2 ). The 
masses of ai and p are easily computed on the same lattice. However, since the range of 
q 2l s that we are probing is very short and since our results are not sufficiently accurate we 
cannot distinguish amongst different interpolating formulas, and we choose to quote our 
value obtained through the linear interpolation to q 2 = 0. All results are listed, in table HI 
and the illustration is provided in fig. [3] for the case of g c (q 2 ,m 2 ). In that same table we 
also give the results for F u ^ d (0, m 2 ) obtained after combining the interpolated values of the 
vector form factors as indicated in eq. 

Second step is far more delicate. In order to reach the physically interesting result 
from our calculation an extrapolation to the chiral limit is required. The relevant formula, 

2 The contribution to the magnetic moment of the heavy-light vector meson that comes from the in- 
finitely heavy quark is zero. 
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p=5.40 




-0.4 -0.2 0.2 

2 2 

q r o 



Figure 3: Linear interpolation in q 2 to q 2 = 0. Illustration is provided with the unquenched data obtained 
at (3 = 5.40 with Nf = 2, which are denoted by the full symbols. The empty symbols show the results of 
interpolation for each of our three light quark (corresponding pion) masses. We stress that the results are 
obtained with the fully unquenched light quark. 



derived in HMChPT, reads [231 El] 

with the particularly interesting parameter g^' which appears in all the HMChPT formulas 
that are being used to guide the chiral extrapolations of the phenomenologically interesting 
quantities in charm physics computed on the lattice. In the expression fl3TT) we multiplied 
with suitable factors of r^/a to make the expression independent on the lattice spacing, a. 
Two important underlying assumptions are made when applying this formula: (i) Eq. (l3~Tj) 
is obtained in the static heavy quark limit and its use is equivalent to the assumption 
that the 1/m™- -corrections do not modify the chiral behavior of this coupling. That 
assumption is quite reasonable because the chiral dynamics is driven by the light quark 
vertex (in our case the current qj^q), while the heavy quark is only a spectator; (ii) 
The HMChPT formula applies to our data too. Of our three light quarks in each set, 
one is heavier than the physical strange quark, one is somewhat lighter and our lightest 
dynamical quark is about a third (half) of the mass of the strange quark when (3 = 5.29 
(5.40). The assumption that such a heavy light quarks obey the behavior driven by the 
HMChPT formula is strong and the check of its validity is one of the most active research 
topics in the lattice QCD community. With these two underlying assumption in mind we 
fit our data to eq. fl3Tj) and obtain 

^HMChPT = {0-63(11)0=5.29, 0.68(7)^ =5 . 4 o} , (32) 



(OK 2 



(4vr/r ) 



'm 7r r ) 2 log(m 7r r ) 2 + 



m. 



(31) 
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while at the physical pion mass we have, g c (m^ hys 2 ) = {0.66(11)^=5.29, 0.71(7) ( g =5 .4o}. If the 
logarithmic correction in eq. (l3Tj) is set to zero we obtain the simple linear extrapolation 
which gives g^ eax = {0. 70(15)^=5. 29, 0. 81(11)^=5. 40 }. As a small check, we combine the 
data of two sets and use only those obtained in the case in which the sea quark mass 
is lighter than the physical strange quark mass, and we obtain fi'HMChpx = 0.67(14), and 
j c (™f ys 2 ) = 0.69(14). These results are to be compared with the quenched value obtained 
in ref. [3], g c = 0.66 (8) (5). We conclude that no effect due to unquenching is visible from 
our results. 

On the other hand, the result for the coupling to the soft photon is completely new. It 
has never been computed even in the quenched approximation. HMChPT formula for this 
coupling was obtained in refs. [2H [25] and reads 



F^m^) = -y/m D m D * 



3" V 2(47r/r )2^ u ' ^ * ^ J Airf 
I f 3 1 + Igj 2 2 \ g\ + d d 



33) 



3^ V 2(47r/r ) 2V ' v ' J 4tt/ 2 

with /3-being the dimension-full parameter which is a measure of the magnetic moment of 
the light quark inside the static heavy-light vector meson, and d u> d are the counter-term 
coefficients which are supposed to be obtained from the fit. As it can be seen from table HI 
our data for Fd(0,m^) are consistent with zero, whereas the results for F u (0,m^) change 
very little when varying the pion mass, contrary to what the HMChPT would suggest. 
In other words, the part linear in m n in eq. (|33|) is far too large to provide a suitable 
description of our data, unless one accepts a huge value for the counter-term coefficient 
{d u ,d) which would then contradict the perturbative feature of HMChPT. For that reason 
we only linearly extrapolate our data and arrive at 

F u (0) e {(5.3 ± 2.3)^ =5 .29, (3.9 ± 1.1)^=5.40} , 

F d (0) e {-(0.1 ± 1.4)^=5.29, -(0.4 ± 0.6)^=5.40} • (34) 

We hope to revisit this issue when more accurate values for F u (m 2 ) obtained with ever 
lighter become available. 



5 Phenomenology 

In this section we discuss our physical results and make a phenomenological discussion 
based on the currently available experimental information. Furthermore, we compare our 
results for the magnetic moment with the results obtained by other theoretical approaches 
such as QCD sum rules and quark models. 



5.1 Our results 

From our results obtained with HMChPT we conclude 

9c G {0.66(11)^=5.29, 0.71(7)^=5.40} g D *Dn+ e {(18.6 ± 3.2)^=5.29, (20.1 ± 2.1)^=5.40} -(35) 
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The expression for the decay width reads 



C 



T(D* -> Dtt) = ^—2- g^nJhf , (36) 

where C — 1 if the outgoing pion is charged, and C = 1/2 if it is neutral. For the charged 
decaying meson we have 3 







) -ml = 


y- \k n -\ 



|jfe | = J( m D* ~ m D + m l 

2m o* 
and we get 

( (72 ± 24) ff - keV, (33 ± 1 1)^ keV , /3 = 5.29 , 
T(D*- ^ Dti) = \ (38) 
[ (82 ± 17) w - keV, (38 ± 8)^ keV , /5 = 5.40 . 

As for the couplings to the soft photon, our results are 

g D *+D+<y = {-0.1(7)^=5.29,-0.2(3)^=5.40} GeV" 1 , 

9Df0D0j = {(2.7 ±1.2)^=5.29, (2.0 ± 0.6)^=5.40} GeV- 1 . (39) 

In other words, the decay width for the radiative decay of a charged D*-meson is negligibly 
small (consistent with zero) whereas for the neutral one we have 



(48±31)keV, = 5.29, 

r(/r° ^ /A) = ^ = <! (40) 

(27 ± 14) keV , (3 = 5.40 , 
where 2m D ,a k 2 , = m 2 D , — m 2 D0 , i.e. fc 7 = 137.1 MeV. 

5.2 Comparison with experiment, other methods and more phe- 
nomenology 

• Let us remind the reader that from the width of the charged .D*-meson measured 
by CLEO, T(D* + ) = 96 ± 22 keV [Tj, and by using the experimentally established 
B(D*+ -> £>+7) = 0.016(4) [25], we get 

r exp '( J D*+) [1-0.016(4)] = T(D* + -> DV) + T(D* + -> D+vr ) 



3 For convenience, in addition to f v = 131 MeV, we list the numerical values of the meson masses used 
in this paper [26] : 

m D + = 1869.6(2) MeV , m D ,+ = 2010.3(2) MeV , m„+ = 139.6 MeV , 
m D o = 1864.8(2) MeV , m D *o = 2007.0(2) MeV , m^o = 135.0 MeV . 
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2m D o\k n + | 3 + m D +\k^o 
127rm D ,+ /2 



-9 2 c 9c = 0.61(7), (41) 



a well known result which is in good agreement with the value obtained in this paper 
on our finer lattice, g c = 0.71(7). 

With the experimental g c and by using the isospin symmetry we can evaluate the 
width of the _D*°-meson as 

T(D*°) = T{D*° -> D-7T+) + T(D*° -> D°tt°) + T(D*° -> L>° 7 ). (42) 

The first channel turns out to be kinematically forbidden. For the last one, instead, 
we have the experimental information that B(D*° — ► -D°7) = 0.381(29) [26J, so that 
with the g c value extracted from eq. (}4"Ij) we have 



1 m D o|A/o| 3 
(0.381 ±0.029) 12vr m D *o * 



r(D* u ) = — — r " ' ^ = 68 ± 17 keV , (43) 



where \k' o\ = 43 MeV. The above result is well below an old -but the only available- 
experimental bound, T(D*°) < 2.1 MeV 



To be able to compare our results for go*D-y couplings with experiment, we first 
combine T(D*+) = 96 ± 22 keV [I], and B(D* + -> D+-f) = 0.016(4) [26], to get 

rv 

B (D*+ -> D + 7 ) x F(D* + ) = -^g 2 D , D+1 k* =► g D , +D+1 = 0.50(8) GeV" 1 , (44) 

where, in this case, fc 7 = 135.6 MeV. Similarly, from the result in eq. (]43p and 
B(D*° -> D°7) = 0.381(29) [26], we obtain 

# D ,o D o 7 = 2.02(26) GeV" 1 . (45) 

We see that these results are in a quite good agreement with values obtained on our 
finer lattice, namely gD*+D+-y = —0.2(3) GeV -1 , and gD* D -y = 2.0(6) GeV" 1 . 

We can also compare to the experimentally measured ratio [261 I2E] 

rm*° _> _d°7t ) , 

i?o = r ^,o ~ D o 7) = L74 ± °- 02 ± °- 13 • (46) 

From our data we have #0 = Llta? at /3 = 5.29, and i? = 2.4±^;| at /5 = 5.40. 
Broadly speaking both our lattice results are consistent with experiment. It will 
be interesting if one can increase the accuracy by which this ratio is computed and 
either check on the validity of various lattice actions that are currently used in the 
literature, or to check on the chiral extrapolations. 
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Method 


Ref. 


T(D*+ -> D + -y) 


r(D*° -> D°-f) 


5-D* + -D+7 




Light Front Model 


m 


0.90(2) 


20.0(3) 


0.4 


1.8 




m 


0.56 


21.7 


0.3 


1.9 


Chiral Quark Model 


m 


1.0 4- 1.5 


38.5 4 43.5 


0.4 4 0.5 


2.542.6 




m 


0.15 4 0.25 


11 4 13 


0.16 4 0.20 


1.44 1.5 


Schrddinger-like Mode7("S") 


m 


0.28 


17.4 


0.2 


1.7 


("V") 




0.6 


14.3 


0.3 


1.5 


Salpeter-like Model 


m 


0.46 


20.8 


0.3 


1.8 


Bag Model 


m 


1.72 


7.18 


0.5 


1.1 


Moment 3 pts -Sum Rule 


m 


0.03(8) 


7.3(2.7) 


0.02(11) 


1.1(2) 


Light Cone Sum Rule (LCSR) 


m 


1.5 


14.4 


0.5 


1.5 


Lattice QCD 


this work 


0.8(7) 


27(14) 


-0.2(3) 


2.0(6) 


Experiment 


m 


1.54(52) 


26(7) 


0.50(8) 


2.02(26) 



Table 5: The table of results for the radiative D* meson decays as computed by various quark models, 
QCD sum rules, lattice QCD and those extracted from experiment in the way discussed in Sec. 5 of this 
paper. "S" /"V" labels the scalar/vector potential in the model of ref. [33]. 



The ratio Rq, as well as the quantities discussed in this section, are interesting because they 
do not involve the weak interactions (CKM matrix elements) and the comparison between 
lattice QCD and experiment requires less assumptions. 

Finally, since our result for the soft photon coupling is new, we should also compare it 
to the existing results in the literature. In table [5] we collect the results obtained by using 
various quark models and various QCD sum rule techniques and compare them to the 
values obtained in this paper. In extracting the couplings gn*D-y from various papers, we 
used eq. (HOI) and the results for the decay widths quoted in the cited papers. To that list 
we should add ref. [38] in which the radiative D*-meson decay was used to fix the unknown 
parameter x(l-3 GeV)<p(l/2), the product of susceptibility of the quark condensate and the 
twist-2 photon distribution amplitude at the middle point. More precisely, the problem of 
large contributions due to coupling to the first radial excitations [39] was circumvented by 
applying the QCD sum rule to the branching fraction, rather than the decay width alone. 

6 Summary and perspectives 

In this paper we report the results of the first lattice QCD study of the gD*D-y coupling, 
i.e. of the magnetic moment of the charmed vector meson. Our result for the neutral D* 
meson is fully consistent with the current experimental decay width T(D* — > D +r y). On 
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the other hand our value for the charged meson suffers from large statistical errors but it 
is small and both consistent with zero and with the measured T(D* — > D +r y). 

We also provide the results of our new computation of the axial coupling go*Dn, which 
corroborates the previous (quenched) result obtained on the lattice. Our value agrees with 
the experimental result for y(D* + ) cleo , with a tendency towards its larger values. 

The results presented in this paper are obtained by using the non-perturbatively 0(a)- 
improved Wilson quarks, both dynamical and valence ones. In spite of the improved 
technique to extract the form factors using the twisted boundary conditions and the con- 
venient double ratios, the gain is not that obvious because our statistical errors are large. 
Besides an obvious goal to reduce the statistical uncertainties one should also try and 
probe the smaller light quark masses and thereby arrive at the precision determination of 
the couplings Qd*d% an d Qd*d^- We argued that the comparison of the accurate lattice re- 
sults for these quantities may be a good testing ground for various lattice QCD approaches 
(different quark actions, effective treatments of heavy quark, chiral extrapolations) against 
experiment, because they involve the calculation of hadronic matrix element but, being the 
strong and electromagnetic processes, no CKM parameter is needed. The approach used 
here could also be used to compute the Qk*k^ coupling for which the experimental value 
is already known to a 5% accuracy. 

Finally the unquenching made here is done with iVf = 2 dynamical quarks and a 
natural step forward would be to include the effects of the 'strange' dynamical quark, 
N { = 2 -> 2 + 1. 
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